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Abstract 



■ The so-called curvaton mechanism -a way to convert isocurvature perturba- 

j_j ■ tions into adiabatic ones- is investigated both analytically and numerically 

in a pre-big bang scenario where the role of the curvaton is played by a suf- 
ficiently massive Kalb-Ramond axion of superstring theory. When combined 
with observations of CMBR anisotropies at large and moderate angular scales, 
the present analysis allows us to constrain quite considerably the parameter 
space of the model: in particular, the initial displacement of the axion from 
the minimum of its potential and the rate of evolution of the compactification 
volume during pre-big bang inflation. The combination of theoretical and 
experimental constraints favours a slightly blue spectrum of scalar perturba- 
tions, and/or a value of the string scale in the vicinity of the SUSY-GUT scale. 



I. INTRODUCTION 



At present, the largest-scale temperature fluctuations of the Cosmic Microwave Back- 
ground radiation (CMBR) are consistent with a (quasi) scale-invariant spectrum of Gaussian 
primordial curvature fluctuations [1-3]. The analysis of the first acoustic oscillations occur- 
ring on shorter angular scales adds the information that such curvature fluctuations should 
be predominantly adiabatic [4-7]. Although sufficiently small amounts of non-Gaussianity 
and/or isocurvature perturbations are not excluded, the above-mentioned observational fea- 
tures represent an important constraint for any scenario trying to model the initial stages of 
our Universe. In a previous Letter [8] we tried to confront the pre-big bang scenario [9-11] 
with these constraints. The present paper contains a full description of that analysis and 
completes it. 

Let us recall that, during the pre-big bang phase, the quantum fluctuations of all the light 
modes present in the low energy effective action are parametrically amplified. Nonetheless, 
sizeable large-scale adiabatic fluctuations are not easily produced from the initial vacuum 
through the usual mechanism of parametric amplification. In particular, both tensor and 
scalar-metric fluctuations are amplified with very steep spectra [12,13], resulting in adiabatic 
modes which are far too small to explain the observed level of large-scale CMBR anisotropics 
[14]. 

However, not all the primordial spectra of pre-big bang cosmology are blue. For instance, 
in a pure gravi-dilaton background, the pseudo-scalar supersymmetric partner of the dila- 
ton in the dimensionally-reduced string effective action, the so-called Kalb-Ramond axion, 
emerges from the pre-big bang phase with a fluctuation spectrum whose tilt depends on the 
rate of change of the compactification volume [15,16]. Depending on this, the axion tilt can 
be negative (red spectrum), positive (blue spectrum), or zero (scale-invariant spectrum). 
However, since the homogeneous (background) component of the axion is trivial, such a 
spectrum does not affect directly the metric, hence no curvature perturbation is generated 
at this primordial level. In other words, the axion perturbations are entropic, isocurvature 
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perturbations. This feature persists, unfortunately, under ^-duality transformations [15,16]. 
If such an axion is massless, or at least light enough not to have decayed yet , the induced 
CMBR fluctuations on large scales can fit the COBE normalization [17-20], but, being nei- 
ther adiabatic nor Gaussian [21], they are not able to fit the observed structure of the first 
few acoustic peaks. 

A possible way out of this problem [8,10,22-24] is offered by the alternative scenario 
of a massive axion, initially displaced from the minimum of its non-perturbative potential. 
In that case axion perturbations couple to scalar metric perturbations through the non- 
vanishing axion's VEV. Eventually, the axion relaxes toward the minimum of the potential 
and then, if heavy enough, decays prior to nucleosynthesis. During the relaxation process 
the dominant source of energy undergoes a drastic change: it consists of the radiation 
produced at the end of the pre-big bang evolution, and later becomes the pressureless fluid 
corresponding to the damped coherent oscillations of the axion. This non-trivial evolution 
results in a non-adiabatic pressure perturbation which, in turn, is well known [25,26] to 
induce curvature perturbations on constant energy (or comoving) hypersurfaces even on 
superhorizon scales. 

The interplay of such different sources of inhomogeneity, throughout the different stages 
of the background evolution, eventually determines the spectral amplitude of scalar curvature 
perturbations right after matter-radiation equality, when all the scales of interest for the 
CMBR data are still outside the horizon. This conversion of isocurvature into adiabatic 
perturbations, originally suggested in a different context by Mollerach [27], also applies to 
more general cases [23,24]. 

Depending upon the initial value o\ of the Kalb-Ramond background, different post- 
big bang histories are possible. If cr < dj C 1 in Planck units (see below, Eq. (3.7), 
for the definition of a), the axion oscillates for a long time before becoming dominant 
and eventually decays. For cr; < a it may never fully dominate the energy density before 
decaying. If, instead, a x ^> 1 the axion will dominate before oscillating and a slow-roll 
(low-scale) inflationary phase could take place in that epoch. As we shall see, of all these 
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possibilities CMBR observations seem to favour the "natural" one, ai ~ 1. In any case, even 
if different post-big bang histories will lead to different spectral amplitudes of the Bardeen 
potential, adiabatic scalar metric perturbations will always be present at some level outside 
the horizon, prior to decoupling. 

The purpose of the present paper is to report on the calculation of the spectral amplitude 
of the induced adiabatic metric perturbations, and on the comparison of the predictions of 
the pre-big bang scenario with the observations coming from the physics of the CMBR 
anisotropies. In order to achieve this goal it is mandatory to have a good understanding 
both of the axion relaxation mechanism and of the evolution of the inhomogeneities. Hence 
analytical results will be supported with numerical examples and vice versa. We will present, 
in particular, a full derivation of the results for the final adiabatic spectrum of the Bardeen 
potential (some of these results have been summarized already in [8]). 

The paper is organized as follows. In Section II the basic equations describing the post-big 
bang evolution of the inhomogeneities and of the background geometry will be introduced. 
In Section III the physics of the different post-big bang histories will be analysed. In Section 
IV the evolution of the background and of its perturbations will be discussed for the case 
in which the amplitude of the initial axion background is smaller than 1 in Planck units, 
Oi < 1. In Section V we will discuss the evolution of the system in the complementary 
case <7i > 1. Section VI is devoted to the phenomenological implications of the large- 
scale adiabatic perturbations produced through the relaxation of the axionic background. 
The obtained results will be compared with observations. Constraints on the pre-big bang 
parameters will be derived. Section VII contains our concluding remarks while, in the 
Appendix, a self-contained derivation of the axionic spectra produced by the pre-big bang 
evolution has been included. 
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II. BACKGROUND AND PERTURBATION EQUATIONS 



As already mentioned, we shall start our analysis at some time 771 in the post-big bang 
epoch, assuming that the axion field has inherited from the preceding epoch appreciable 
large-scale fluctuations, while other sources of energy as well as the metric are exactly 
homogeneous. It will also be assumed that, initially, the dominant source of energy is in the 
form of radiation. The post-big bang dynamics takes place, in the present analysis, when the 
curvature scale has fallen to a sufficiently small value (in string units) so that the use of the 
low-energy effective action is appropriate. Furthermore, for 77 > rji the dilaton is assumed to 
be frozen already at its present value. 

Under these assumptions, the evolution of the geometry is determined by the Einstein 
equations, supplemented by the conservation equations determining the dynamics of the 
sources*: 

K-\€R = \(n{a) + T£), (2.1) 
8V 

<7 a/3 V Q V /3 <7 + — = 0, (2.2) 

where T^(a) and 7~f are, respectively, the energy- momentum tensors of the axionic back- 
ground and of the matter fluid. Notice that the covariant conservation of T^(cx) is dynam- 
ically equivalent to the evolution equation of the axionic field, i.e. Eq. (2.2), and implies, 
through the contracted Bianchi identities: 

v Q r/ = 0. (2.3) 

In a conformally flat background geometry, 

ds 2 = a 2 (r])[dr] 2 - dx% (2.4) 



* Gravitational units 16ttG = 1 will be used throughout. When explicitely written in the formulae, 
Mp = (167rG) _1 / 2 = 1.72 x 10 18 GeV. In these units, a is the canonically normalized axion field. 
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Eqs. (2.1)-(2.3) lead to a set of three independent equations, whose specific form is dictated 
by the fluid content of the primordial plasma. In the case of a radiation fluid we have 

%° = Pr, V = -p T 5l p r = |, (2.5) 

and Eqs. (2.1)-(2.3) lead to 
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a' 2 

Pr + — ~V 

a 1 



(2.6) 



a" + 2Ha> + a 2 ^- = 0, (2.7) 
da 

p' x + AHp r = 0. (2.8) 

Here the prime denotes the derivation with respect to the conformal time coordinate 77, and 
H = (In a)'. For future convenience we also recall that the connection between Ti and the 
Hubble parameter is H = H/a. The effective energy and pressure densities of a will be 
given by 

P. = £ + V, V. = £ - V. (2.9) 



The set of dynamical equations (2.6)-(2.8) is supplemented by the Hamiltonian constraint 

a 2 r a' 2 1 

which imposes a specific relation on the set of initial data and is required, in particular, for 
the numerical integration of the background evolution. 

During the post big-bang phase, the first order perturbation of Eqs. (2.1)-(2.3) provides 
the linear (coupled) system of evolution equations of the inhomogeneities. To first order in 
the scalar metric fluctuations, the line element (2.4) can be written as [26] 

ds 2 = a 2 (r]){(l + 2(j))drf - 2d i Bdx i dr ] 2 - [(1 - 2^)8^ + 2d i d j E}dx i dx j }. (2.11) 

Since there are two gauge transformations preserving the scalar nature of the above metric 
fluctuations (ip,<j>,E,B), two gauge-invariant (Bardeen) potentials can be defined [26,28]: 
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$ = <j> + -[(B-E , )a]', 
* = 1 f> -H(B - E'). 



(2.12) 
(2.13) 



Appropriate gauge-invariant variables can also be defined for the perturbations of the 
sources, in such a way that 



x 0p) = 5a + a'(B-E'), 
5pW> = 5p I + f/ I (B-E'), 



,(gi) 



vr'=v t + {B-E'), 



(2.14) 
(2.15) 
(2.16) 



whose physical interpretation is particularly simple in the so-called longitudinal gauge [26] in 
which E = = B. Here 5Tq = 5p r) and the velocity potential is defined by the off-diagonal 
fluctuations of the radiation energy-momentum tensor as 



5%° = (p r + p x )v?5u h 



(2.17) 



where u° = 1/a and, in the longitudinal gauge, 5ui = adiV r . 

By perturbing the diagonal components of T"(cr), and using Eqs. (2.12) and (2.14), 
the fluctuations of the axionic energy and pressure densities can be expressed in a fully 
gauge-invariant way as follows t; 



Spa = — 

a 1 

5 Pc = "5 



-<5>a' 2 + a' X ' + -7^a 2 X 
-$a' 2 + a' X ' - —a\ 



(2.18) 
(2.19) 



The variables characterizing the gauge-invariant fluctuations of the sources can be defined 
in different, but equivalent, ways [25,29]. For instance, it is sometimes useful (especially 
in the case of fluids with constant speed of sound) to write equations for the combination 
(Sp r /p r — 4$), whose evolution greatly simplifies at large scales. 



tin the following, since we will be dealing only with gauge-invariant quantities, the superscript 
'(gi)" can be consistently dropped without confusion. 
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The fluctuations of the off-diagonal (space- like) components of Eq. (2.1) imply that 
$ = Hence, in terms of the variables defined in Eqs. (2.12)-(2.19), the (00) and (Oi) 
components of the perturbed Einstein equations (acting as Hamiltonian and momentum 
constraints for the evolution of the Bardeen potential) can be written in terms of the gauge- 
invariant velocity potentials v T , v a , and of the radiation and axion density contrasts 5 T = 

Spjp r , 5 a = 5 Pa/ Pa, aS follows: 

,2 



V 2 $ - 3H(H$ + = j (pA + PaSa) , (2.20) 

(p r +p r )v r + (pa +Pa)Va , (2.21) 



a 2 

H<$> + $' = — 
4 



Here the axion velocity potential, v a , is defined by 



x (2.22) 



a^/Pa + Pa 

and is the axionic counterpart of the velocity potential introduced for the radiation fluid. 

The constraints (2.20) and (2.21) are to be supplemented by the dynamical equations 
coming from the perturbation of the (ii) components of Einstein's equations (2.1), of the 
axion equation (2.2) and of the continuity equation (2.3). For the gauge-invariant quantities 
defined above, such dynamical equations are, respectively: 

2 2 

$" + m& + (H 2 + 2H')$ = ^-pA + ^-5pa, (2.23) 

X " + 2H X ' - V\ + ^a 2 X - 4a'& + 2^a 2 $ = 0, (2.24) 

b[ - 4$' - \v 2 v T = 0, (2.25) 

v' t - ^5 r - $ = 0. (2.26) 

Finally, the perturbation of the covariant conservation of the axionic energy-momentum 
tensor leads to two useful equations: 

pX - (Pa + Pa) VV - 3H P a5a - 3$ ' { Pa + Pa ) + 3H5 Pt7 = 0, (2.27) 

v'a + Un + P ^Av a - - $ = 0, (2.28) 

V Pa+Pa-> Pa + Pa 
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which are implied, as it should be, by Eqs. (2.23)-(2.26) when the background equations 
(2.6)-(2.7) are used. 

It is also useful to notice that, by combining Eqs. (2.10) and (2.23), we can eliminate 
the fluid variables, and we obtain 



which, together with Eq. (2.24), provides a closed system of equations for $ and \- Of 
course, the velocity potential and the density contrast of the fluid do not disappear from 
the physics of our problem, and have to be directly computed using the Hamiltonian and 
momentum constraints of Eqs. (2.20) and (2.21). 



Given the system of Eqs. (2.23)-(2.26), supplemented by the constraints (2.20)-(2.21), 
it is sometimes appropriate to select variables obeying simple evolution equations in the 
long-wavelength limit, in which the spatial gradients are negligible. For this purpose, a 
particular combination of Eqs. (2.20) and (2.23) will be considered, and the fluctuations in 
the total energy and pressure densities will be defined: 



In terms of the quantities defined in Eq. (2.30), the evolution of the Bardeen potential can 
be formally written in terms of a single equation 




(2.29) 



A. Curvature perturbations from non adiabaticity 



(2.30) 



$" + 3H{1 + c£)$' + [2H' + H 2 {1 + 3c")]$ - c^V 2 $ = ^[Sptot - c 2 s 5p tot ], (2.31) 



where c s is the speed of sound for the total system, defined by 



c 2 = Kot = P'a+Pr 
S P'tot Pa + Pr 



(2.32) 



or, using the explicit form of the background equations, 
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.2 



c: = + + (2.33) 

3 I Pr + f(p«x + Pa) J 



where V^ CT = dV/da. 

The left-hand side of Eq. (2.31) (except for the Laplacian term) can now be expressed 
as the time derivative of a single gauge-invariant function (, namely 



c = - 



. + f*(»*±*)] S -(. + H™±£), (,34) 
a 2 vptot + Ptot / J ' H 1 -W J 



where the second equality follows by using the background equations of motion (2.6)-(2.10). 
By using this variable, Eq. (2.31) can be written as 



d( H . mc 



,2 

^Pnad " 2( "7 , V 2 $. (2.35) 

dr] p tot + p to t a 2 (ptot + Ptot) 

where we have defined 

^Pnad = Sptat ~ C^p t ot- (2.36) 

As noticed long ago [25,28,30], the variable ( represents the inhomogeneities in the spatial 
part of the space-time curvature, measured with respect to comoving hypersurfaces (a = 
constant). Using Eq. (2.21), the variable ( can also be usefully related to the total velocity 
potential as 

C=-($ + Wutot), (2-37) 

where 

(Ptot + Ptot Mot = (Pr + PrK + (Pa + P a )v a . (2.38) 

In our specific case, using the full set of background and perturbation equations in the long- 
wavelength limit, where 5 T ~ 4$ according to Eq. (2.25), the expression for ( can be written 
in the following convenient form 

C " 4p r + 3(p CT + Pa ) ■ ^ 



As we will discuss in detail in Sect. V, in the absence of a dominant radiation fluid 
<5Pnad is zero at large scales, i.e. up to terms containing the Laplacian of <£>. However, in 
a radiation dominated regime, <5p na d 7^ and Eq. (2.35) implies (' 7^ even in the long- 
wavelength limit. Let us then compute the general form of 5p na d, for the full system of axion 
plus fluid perturbations. By using the previous definitions we obtain 

^nad =Pr(l~ 4) *r + *{£ ~ l)(p. + Pa) + ^(1 " # - ^ X (l + (2.40) 



On the other hand, using Eqs. (2.20) and (2.21), we can write 

4 u'y rcrV dV 

®(P* + Pa) = "-V 2 $ + m-±- + AHfhVr + pA + ~t + JT X ■ ( 2 - 41 ) 

(X (X L CL CJCT 

Thus (neglecting the spatial gradient of $) we get 



<5pnad = -IpA - + 4H P Mc 2 s - 1) + 3ft(^ - 1)^. (2.42) 

d OCT GT 

The above equations are useful to compute, in some specific phase of the dynamical evolu- 
tion, the source term of Eq. (2.35) whose integration allows us to obtain the explicit time 
dependence of (. In [8] we have determined the evolution of the fluctuations by following 
the ( variable. In the present investigation we will solve the perturbation equations both in 
terms of <3> and (, checking numerically the consistency of the two approaches. 



III. POST BIG-BANG HISTORIES 

At the beginning of the post-big bang evolution the background is characterized by a 
"maximal" curvature scale H 1 , whose finite value regularizes the big bang singularity of the 
standard cosmological scenario, and provides a natural cutoff for the spectrum of quantum 
fluctuations amplified by the phase of pre-big bang inflation (see below, in particular Section 
VI) . In string cosmology models such an initial curvature scale is at most of the order of the 
string mass scale, i.e. H 1 < M s ~ 10 17 GeV. 

The Kalb-Ramond axion has gravitational coupling to photons and to the QCD topolog- 
ical current but it is not necessarily identified with the invisible axion [32] usually invoked in 
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the explanation of the strong CP problem via an initial misalignment of the QCD vacuum 
angle i? [33]. The potential of Kalb-Ramond axion is, strictly speaking, periodic. The peri- 
odicity of the potential occurs whenever a Peccei-Quinn symmetry is spontaneously broken 
down to a discrete symmetry corresponding to shifting the i? angle by multiples of 2n (see, 
for instance, [34]). However, close to the minimum of the potential (i.e. sufficiently late in 
the process of relaxation) the potential can be assumed to be quadratic. Such an approxi- 
mation is expected to be realistic for values of a that are small compared to its periodicity. 
Unfortunately, translating periodicity in d into periodicity in a involves a normalization 
factor that is unknown in the strong-coupling region where the dilaton is supposed to be 
frozen at late times. For this reason, we shall keep the initial dispacement in Planck units, 
o"i, as a free parameter. 

We start our study of the background and perturbation evolution at an initial curvature 
scale Hi < Hi, when the energy density of the background is mainly stored in the radiation 
fluid, while the energy density of the axion is dominated by the potential: 

Pl (vi)^>Pa(m)-V(vi). (3.1) 

During the first stages of the evolution a remains approximately fixed at the initial value a- x 
up to corrections 0(V a ). In the course of such a "slow-roll" phase, the curvature scale of 
the background decreases, until it becomes comparable with the curvature of the potential. 
The axion background will then start oscillating, at a typical scale 

H osc ~ m (3.2) 

(note that, as already mentioned, we are assuming that the potential is quadratic). At the 
curvature scale 

H a ~ma(t), (3.3) 

the axion field will dominate the background. The specific value of the scale H a depends 
upon <7i and also upon the evolution after r\i. In fact, during the oscillatory phase the axionic 

11 



energy density decreases, on the average, as a~ 3 , i.e. slower than the energy of the radiation 
background, p r ~ a -4 . From Eqs. (3.2) and (3.3) it is then clear that, depending on the 
initial value of a, the oscillations of the axionic background may arise either before or after 
the phase of cr-dominance. 

Irrespectively of its initial value, the coupling of a to photons is gravitational, i.e. sup- 
pressed by the Planck mass. The decay takes place when the curvature scale is of the same 
order as the decay rate, namely when 

Try-* 

The late decay of a is in general associated to a significant entropy release, which has to 
be carefully constrained [35-37] not to spoil the light nuclei abundances and the baryon 
asymmetry generated, respectively, by primordial nucleosynthesis and baryogenesis. 

In our context, for typical values of Hi, and for a realistic scenario, the decay of a is 
constrained to occur prior to nucleosynthesis, i.e. at a scale H d > if N ~ (lMeV) 2 /M P , which 
implies m > 10 TeV. The lower bound on the axion mass is even larger if we require that the 
decay occurs prior to baryogenesis at the electroweak scale (characterized by a temperature 
of the cosmological plasma of order 0.1 TeV), which implies m > 10 4 TeV. If, on the contrary, 
baryogenesis occurs at a large enough scale preceeding the phase of axion dominance and 
decay, then the minimal value of m allowed by the entropy constraints [35-37] is, in general, 
Ui-dependent. In that case, however, the resulting lower bound is strongly dependent on the 
given model of baryogenesis, and can be somewhat relaxed by various mechanisms. In the 
rest of this paper we will thus adopt a conservative approach, by taking the nucleosynthesis 
bound m > 10~ 14 M P as a typical reference value. 

A. Late dominance of the axion: o\ < 1 

If a(r)i) — <Ji < 1, then the axionic background first experiences a phase of radiation- 
dominated oscillations, from H osc down to H a . The duration of this phase depends upon 
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<7i, since (a osc /a a ) ~ of, and it may be rather long, if o- x <C 1. During this phase the axion 
potential energy decreases as a~ 3 . Consequently, the typical scale of axion dominance is, 
from Eq. (3.3), 

H (j ~ ma?. (3.5) 

From H a down to H d: i.e. inside the regime of axion-dominated oscillations, the effective 
equation of state of the gravitational sources, averaged over one oscillation, mimics that of 
dusty matter, with (p a ) = 0. During this regime the scale factor and the Hubble parameter 
also have oscillating corrections, which vanish on the average, and decay away for large times. 
It should be stressed, however, that the effective equation of state of the axion background, 
for curvature scales smaller than H a , depends upon the curvature of the potential around 
the minimum. If, for instance, the potential is not quadratic, but quartic, the coherent 
oscillations will lead to an effective equation of state that simulates a radiation fluid, i.e. 
3(p CT > = (Pa) [38]. 

The occurrence of the axion-dominated phase requires 

H d < H a , (3.6) 

which imposes a lower bound on the initial axionic amplitude, namely 

1 xn > y/m/Mp = a. (3.7) 

This constraint, however, is not so demanding, given the generous lower bound on m (in 
Planck units) allowed by nucleosynthesis and baryogenesis. Finally, after the axion decay, 
the Universe enters a subsequent radiation-dominated epoch. From this moment on, the 
evolution of the background fields is standard. 

B. Early dominance of the axion: a\ > 1 

If <r(?7i) = o"i > 1, then the axion, right after the onset of the radiation-dominated epoch, 
starts again rolling down its potential. This initial part of the evolution is completely 
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analogous to that of the a\ < 1 case. However, for o\ > 1, the axion dominance will occur 
before the onset of the axion oscillating phase, i.e. 



H osc < H a , (3.8) 



where, for a generic potential, 



H a ~ v^R) (3-9) 



(since the kinetic energy of the axion is negligible during the slow-roll evolution). At H = H a 
the Universe enters a phase of accelerated expansion (slow-roll inflation) whose duration, 
for a quadratic potential, is given by 



ry a fina i 

Zj a = = exp 

^initial 



g (^initial °final) 



(3.10) 



This inflationary phase will last until H = H osc ~ m, <Tfi na i ~ 1 (if we assume, again, that 
close to its minimum the potential is quadratic). For H < H osc the background will be 
dominated by the coherent oscillations of the axion, whose decay will eventually produce a 
second radiation-dominated phase (in full analogy with the case a- x < 1). 
This scenario requires, for consistency, that 

H a < Hi < H 1 < M s , (3.11) 

which, in the case of a quadratic potential, amounts to require 

H 1 /m>H i /m>a i >l. (3.12) 

As in the case of Eq. (3.7), this bound is not so restrictive, given the limits on the axion 
mass. Indeed, in the case <7\ > 1, the most stringent constraints are coming not from 
the evolution of the background geometry but, as we shall see, from the evolution of the 
fluctuations that forbid too large values of ov One is then left with a situation where a x ~ 1 
and H a ~ H osc . In such a case, the phase of axion-dominated oscillations will take place 
right after the radiation-dominated period of slow-roll, without a long intermediate epoch 
of inflation. 
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C. Initial conditions for the fluctuations 



Given the coupled system of gauge-invariant perturbation equations, the initial condi- 
tions for the Bardeen potential, for the perturbed radiation density and for the radiation 
velocity field, will be imposed as follows 

$kfai) = 0, 5 r (r)i, k) — 0, v I (r H ,k) = 0, (3.13) 

assuming that no appreciable amount of adiabatic metric perturbations has been directly 
generated (on large scales) by the pre-big bang dynamics. The only non-vanishing initial 
fluctuations are the (isocurvature) axionic seeds, amplified from the vacuum during the 
pre-big bang evolution: 

Xk(m) = Xi(k) ^ 0, (3.14) 

so that, from Eq. (2.18), 5 a (r)i,x) ^ 0. 

In the present analysis we shall assume that the amplitude of the initial axion fluctuations 
is smaller (for all modes) than o - ;, i.e. 

k 3 / 2 \ Xi (k)\ < a,. (3.15) 

In the opposite case, k 3 ^ 2 \xi(k)\ > <7i, we are led to the case already analysed in [17-19,21] 
where a x was assumed to be zero, and the obtained large-scale fluctuations are of the isocur- 
vature type, and strongly non-Gaussian. If a x > a, the non-Gaussianity is rather small, but 
can be enhanced if the axion does not dominate at decay (a - ; < a) [23,39]. 

IV. BACKGROUND AND PERTURBATION EVOLUTION FOR a, < 1 

In view of the forthcoming phenomenological applications, the main quantity that we 
need to compute is the spectral amplitude of the Bardeen potential after axion decay, during 
the subsequent radiation-dominated evolution, as a function of the spectral amplitude of the 
axion fluctuations amplified by the phase of pre-big bang inflation. It is important, for this 

15 



purpose, to have a reasonably accurate control on the evolution of the background and 
of the fluctuations. Using different approximations, motivated by the hierarchy of scales 
discussed in the previous section, we will first analytically determine the evolution of the 
system through the different cosmological stages. Numerical integrations will then be used 
in order to check the analytical results in the cross-over regimes connecting the different 
phases of the background evolution. 



A. The radiation-dominated slow-roll regime 

During the first stage of evolution, p r (r]i) 3> V(r]i). In this limit, Eqs. (2.6)-(2.8) and 
(2.10) simplify to: 

n 2 = j Pl , (4.1) 

p' I + m Pl = 0, (4.2) 

a" + 2Ha' + a 2 ^- = 0. (4.3) 
do 

Equations (4.1) and (4.2) imply that Ha is constant. Furthermore, since the kinetic energy 
of a is subleading with respect to the potential, the axionic field slowly rolls down the 
potential. In such a situation a systematic expansion in the gradient of the potential, V a , 
can be developed, and the background evolution is adequately described by the following 
approximate equation 

✓ = -\&, (4-4) 



which can be integrated to give 



rr-LfYtL-Yf. 
1 20V# 2 H 2 



), (4-5) 



i.e. a is approximately constant up to corrections that depend upon the specific form of the 
potential, and which induce a slight decrease of the axion background. 



16 



In order to solve the Hamiltonian constraint (2.20) it is now convenient to work in terms 
of the Fourier components of the perturbation variables, $k,5 r (k), and so on. Since we are 
interested in large scale inhomogeneities we first obtain, from Eq. (2.25), 



S T (k) ~ 4$ fc , 



(4.6) 



where the integration constants vanish because of Eq. (3.13). Consequently, using Eq. 
(2.18), the Hamiltonian constraint (2.20) can be written as 



-m(n<$> k + $' fc ) - $ k 



/2 

2 U 


1 " 


a p r - — 


" 4 . 



v'x'k + V a a 2 Xk 



(4.7) 



where the spatial gradients have been consistently neglected. Using Eq. (4.1), 

1 



vm 



&Xk + -Q^a Xk 



(4.8) 



On the other hand, from Eq. (2.24), the evolution of Xk is approximately given by 

' ~ - l -v a ! 

Xk — g * / ,cra y_^Xh- 



(4.9) 



The first term on the r.h.s. of Eq. (4.8) thus contains three derivatives of the potential, and 
it is subleading with respect to the second term. Direct integration of Eq. (4.8) then gives 

1 



VrXk + 0(Vi). 



(4.10) 



(4.11) 



As a consequence, from Eqs. (2.25) and (2.26) we can determine <5 r and v x as 

Inserting now the obtained solutions in the remaining equations (2.21) and (2.23) (not used 
for the above derivation), we can see that they are satisfied with the same accuracy. 

The time evolution of Cfc i n the radiation-dominated, slow-roll regime can finally be 
determined from Eq. (2.34): 



(4.12) 
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so that $fc and (k obey the following simple relation 

Mv)^-(V7)Ck(v) + 0(V 2 a ). (4.13) 

It should be appreciated that Eq. (4.12) can also be obtained by direct integration of 
Eq. (2.35). In the limit (p a + p a ) <C p r , Eq. (2.33) implies indeed 

c H + 2^ v - < 414 » 

On the other hand, from Eqs. (2.40) and (4.6), the approximate form of 5p na d(/c) is 

5 Pnad (k)~-^V aXk + 0(V*) (4.15) 

(again, terms with more than one derivative of the potential have been neglected). By 
inserting this result into Eq. (2.35) we are led to the equation 

= (4.16) 
a ma p r 

whose direct integration (recall that p r ~ a -4 ) leads, as expected, to Eq. (4.12). 

The above approximate results for $ fc and ( k hold for a generic (flat enough) potential. 
However, in order to check the correctness of our approximations numerically, it is useful to 
consider the simple case of a quadratic potential: 

V{a) = (4.17) 

In that case, Eqs. (4.4), (4.9), (4.10), (4.12) lead to 



a (t ~ <7; 



i - ^ (- 4 - 1) + 



Xk{r) ~ Xi(fc) 



£ (r 4 - 1) + 0( A 8 ) 



(4.18) 
(4.19) 



Mr) - [S (r 4 - 1) + 0(A 8 )] , (4.20) 

COO - ^ [p 2 (r 4 - 1) + 0(A 8 )] , (4.21) 

where X\(k) is the initial spectrum of the axionic fluctuations, and we have defined the 
(dimensionless) rescaled mass and conformal time coordinate: 
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1 



/i — mrjidi — m/Hi. 



(4.22) 



The time rji is the initial integration time and a; the initial normalization of the scale factor. 
These rescalings are useful in order to compare the numerical results with the analytical 
calculations. 

We have performed a numerical integration by choosing initial conditions at sub- 
Planckian curvature scales, i.e. 



(in Planck units), and setting a\ = 1. Given a value of a x compatible, for a given mass, with 
Eq. (3.7), the constraint (2.10) fixes the initial radiation background p x {j]\) to a value much 
larger than the axionic potential. Similarly, the initial values of the derivatives of and Xk 
are obtained by imposing, on the initial data (3.13)-(3.15), the Hamiltonian and momentum 
constraints of Eqs. (2.20) and (2.21). It has been checked that all the constraints (both for 
the background and for the fluctuations) are satisfied at every time all along the numerical 
integration. The system describing the evolution of the fluctuations, in particular, can be 
integrated in two different (and complementary) ways. We could either use Eqs. (2.23)- 
(2.26) and follow the evolution of all variables, or use Eqs. (2.24)-(2.29) and integrate the 
system only in terms of and Xk- We have performed the numerical integration in both 
ways, and checked that the results are the same. 

In Figs. 1, 2 and 3 we report, as full curves, the results of the numerical integration for 
a quadratic potential. The analytical results of Eqs. (4.18)-(4.21), based on the slow-roll 
approximation, are also illustrated, for comparison, by the dashed curves. 



During the radiation-dominated regime, and for a quadratic potential, the axion evolu- 
tion equation (2.7) can be written as 



Hi 




m > i 



(4.23) 



B. Radiation-dominated oscillations 




g = aa, 



(4.24) 
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FIG. 1. The full curve shows the result of a numerical integration for the case Hi = 0.01, and 
for the set of parameters reported in the figure. The dashed curve shows the approximate analytical 
solution based on Eqs. (4.10) and (4.20). 
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FIG. 2. Evolution of Xk, reported for the same set of parameters as in Fig. 1. The dashed 
curve corresponds to the approximate analytical result of Eq. (4.19). 
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FIG. 3. Evolution of Cfc, for the same set of parameters as in Figs. 1 and 2. The numerical 
result (full curve) is compared with the approximate analytical result of Eq. (4.21) (dashed curve). 



and its exact solution can be given in terms of Bessel functions [40] as 



[XT 



(4.25) 



where C1/4 is a linear combination (with constant coefficients) of Bessel functions of order 
1/4 and fir 2 /2 ~ m(t-ti). By imposing the correct boundary conditions and normalization, 
in such a way that er(r) — > a x for r — > 1, we obtain 

1 



(TIT) 



0\ 



-Jl 



(V). 



Vi < V < Va, 



(4.26) 



where J1/4 is the first-kind Bessel function. Notice that the small argument limit of this 
equation, for /j, 1 and r — > 1, exactly gives the result (4.18), obtained in the slow- roll 
approximation. This exact analytical solution can also be used as a consistency check of 
the quadratic approximations when the potential, during slow-roll, has a more complicated 
analytical form. 

The onset of the axion oscillations can be determined from the first zero of J\/a([it 2 /2), 
which occurs for 



2.78, 
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(4.27) 



namely for 

r m = ^, ei ~ 2.35. (4.28) 

Different definitions of the oscillation starting time, for instance related to the breakdown 
of the slow-roll approximation, would lead to similar numerical values, i.e. to Eq. (4.28) 
with ei ~ (12) 1 / 4 . In the large argument limit (/i/2r 2 1) of Eq. (4.26) the solution finally 
describes the oscillating regime, 



where the phase and amplitude of oscillations are fixed by the initial conditions. 

An approximate solution of Eq. (2.24), similar to Eq. (4.26), holds for the axion fluctu- 
ations, namely 

Xk{r)^ Xi{ ^ - ) — J^fgr 2 ), Vi<V<V*, (4.30) 



V? J 1/4 (/i/2) 1/4 V 2 
whose small and large argument limits lead, respectively, to 



Xk(r) - Xi(k) 



1 - ^ (r 4 - 1) 
20 V ; 



(4.31) 



We notice that Eq. (4.30) is obtained by solving the approximate evolution equation 

xl + 2H X ' k + m 2 a 2 Xk ~0, (4.33) 

i.e. neglecting the terms containing the Bardeen potentials in Eq. (2.24). In the slow-roll 
approximation, as previously stressed, these terms can be neglected for a generic potential 
term. However, they can also be neglected in the oscillating phase, provided the potential is 
well approximated by a quadratic form. We have explicitly checked that the exact analytical 
solutions (4.26) and (4.30) are in perfect agreement with the results of a numerical integration 
performed with a quadratic potential. 
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FIG. 4. The exact, numerical evolution of a and Xk (full curves) is compared with the inter- 
polating solution (dashed curves) obtained by matching the slow-roll solutions (4.18) and (4.19) 
with the WKB approximated solutions (4.34) and (4.35), valid during the radiation-dominated 
oscillations of the axion. 

Thus, for a potential which is generic during the slow-roll phase (but still quadratic 
during the oscillating regime) it will be sufficient to work out the slow-roll solutions specific 
to that potential, from Eqs. (4.4), (4.9), and match them (with their first derivative) to the 
WKB solutions of Eqs. (4.24) and (4.33), namely 



The matching will allow a determination of the precise amplitudes and phases of a (and Xk) 
in terms of a; (and Xi(^))- 

As an application of this technique let us consider the example of the quadratic potential, 
using the slow-roll solutions for r < r m . The result of this exercise is reported in Fig. 4 
where, with the full curves, we illustrate the numerical results (coinciding exactly with the 
analytical solutions). With the dashed curves we show the interpolating solutions obtained 




(4.35) 



(4.34) 
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by matching Eqs. (4.18) and (4.19) (obtained in the slow- roll approximation) with the WKB 
solutions (4.34) and (4.35), valid in the oscillating regime. 

The time evolution of a{rj) and Xk(v) explains why, for r > r m (i.e. after the slow- roll 
regime where $fc ~ a 4 according to Eq. (4.10)), the Bardeen potential enters a phase of 
linear evolution (in conformal time). This feature is illustrated in Fig. 5, where we report the 
numerical results for the evolution of the Bardeen potential, computed for different values 
of the axion mass. 

An analytical estimate of the slope of the linear regression for after the end of 
the radiation-dominated slow roll, can be obtained from the Hamiltonian constraint (2.20), 
which can be recast in the following form: 



<9 r 4 

1 ' • ; (4.36) 



(r 3 $ fc 



dr K ^ 12 



dcrdxk , 2 2 

+ fx t axk 



r^ 2/2 dx i 

-12 cos (27 + fir 2 ) + 18 / —cos (a; + 7) , (4.37) 

J x J 



.dr dr 

where we only assumed a quadratic form for the axion potential. By using the WKB solutions 
(4.34) and (4.35) we obtain 

9Qaf ^ 384af/ir 3 

where the integral can be expressed in terms of Ci(iu) = — £3 ^ £ dx and Si(iu) = J™ ^f-dx, 
and we have assumed (5 — 7. The oscillating terms are suppressed by r~ 3 and can be 
neglected (in agreement with the numerical results of Fig. 5), since we are considering the 
regime r > r m ^> 1. On top of the oscillating terms, the amplitude of the term responsible for 
the linear growth can be extracted from the numerical solutions by fitting their asymptotic 
behaviour with the line 

$k(v) = o"iXi[-e2 - e 3v ^Zr], r > r m . (4.38) 

In the case of a quadratic potential we obtain 

e 2 ~ 0.001, e 3 ~ 0.0437. (4.39) 

The accuracy of this result can be appreciated from Fig. 5, where the dashed lines (barely 
distinguishable from the numerical solutions) are plotted according to Eqs. (4.38), (4.39). A 
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FIG. 5. The results of the numerical integrations for the Bardeen potential are illustrated by 
the full curves for a quadratic potential and for different values of the mass. The dashed lines 
represent the linear fit of Eq. (4.38). 

different form of the potential will not affect the angular coefficient of the regression (which 
is determined by the phase of the radiation-dominated oscillations), but only the constant 

£2- 

It may be interesting to look also at the analytical estimate of €3, for a quadratic potential. 
The initial amplitudes <r 2 and X2(k) of Eq. (4.37) are determined by the large argument 
limit of the exact solutions, Eqs. (4.29) and (4.32). In this case we get 



in excellent agreement with (4.39). For a generic potential, we could also determine €3 by 
adopting an approximate procedure, i.e. by taking the slow-roll solutions for Xk and a from 
Eqs. (4.4) and (4.9), and matching them in r m to the WKB solutions (4.34) and (4.35), in 
order to determine amplitude and phase. 

The linear growth of the Bardeen potential continues until p a , decreasing as a -3 , equals 
p T . This happens at a time r a such that 




(5/4) 



0.0435, 



(4.40) 
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12H Z ~ a'mV. (4.41) 
The expansion of Eq. (4.29) for \x <C 1 then gives 

From Eq. (4.38) we can then finally obtain the value of the Bardeen potential, at the onset 
of the phase of cr-dominated oscillations. By using the value of r CT given by Eq. (4.42), the 
result is 

<S> k { Va ) ~ e 5 ^, e 5 = e 4 e 3 ^ 0.25. (4.43) 



C. The axion-dominated oscillations 

Using standard techniques suitable for the oscillating regime [38,26], Eqs. (2.10) and 
(2.7) can be solved and, in the case of a quadratic potential, the oscillating terms lead to 
a geometry that reproduces (but only on the average) a matter-dominated Universe. The 
oscillating corrections will be suppressed for large (cosmic or conformal) times, and can be 
easily computed in the cosmic-time gauge, where: 

1 



H 2 



1 



• 2 i 2 2 

a + m a 
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a + 3H& + m 2 a = 0. (4.44) 

Using the auxiliary variable Z = da/dlna, which satisfies 

dZ ( Z 2 \f Z\ . 

A? = - 3 ( 1 - 12)( 1 + 4 ^)' (445) 

and defining two angular variables (r,6), 

Z = VV2 sin 6, a = r cos 6, (4.46) 



(such that H = mr/y/VZ), the following two equations are obtained: 
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\/3 \/3 
9 = — — mr sin 29 — m, r — — —mr 2 sin 2 9. (4.47) 

They can be solved, and the solutions expanded for large times at any order in 1/t. 

A similar procedure can be carried out in conformal time. Equations (2.7)-(2.10) are 
equivalent to the following set of equations 



r' = - 



^-marfr 2 sin 2 ^ ma f ri j ? a = r cos ^ ma ^ r] j (4.48) 



(where a/ is an appropriate dimensionful integration constant), and their solution leads to 
the expansion 

3 (2 maf7] 3 \ 1 



/ \ r 2 ( * TTUlfV \ ,nl 



5- 



2a\m ln if \ 3 / 77 

«(„M? + ^ sin (^!) + 0( i)l, 

Ir] ma^rf V 3 / 77 J 
aM ^^L 3Cos r^f) +0( l). (4.49) 

A posteriori, as a cross-check, Eqs. (4.49) can be inserted into Eqs. (2.6)-(2.10), to see that 
all terms up to 0(1/ if ) cancel, as expected. 

In the phase dominated by the oscillating axion the effective gravitational source is 
pressureless, on the average. By inserting the condition (p a ) = into the background and 
perturbation equations, we get 

(5 a (k)) ~ -2(<t> fe ), (5 r (k)) ~ -2(5 a (k)), (4.50) 

which can be inserted into Eq. (2.39), obtaining 

(Ck)*l(8*(k))~-l(Q k ). (4.51) 

when (p a ) ~ and p a ~^> p r . 

Since p a vanishes only on the average, more accurate solutions have to be supplemented 
by oscillating corrections. Using Eqs. (4.49) an approximate form of the perturbations in 
the oscillating regime can be obtained. One finds that 5 a (k) and are almost constant (up 
to oscillations), i.e. 
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where 



*o(*0 = <**fa)>, 



V* <V < Vd, 



and 



(4.52) 
(4.53) 

(4.54) 



(4.55) 



XoW = -^ $ o(A:). 

The above solutions satisfy the evolution equations of the fluctuations up to 0(r]~ 5 ). 



(4.56) 



D. The axion decay and the subsequent radiation-dominated phase 

When the decay rate of the axion equals the cosmological expansion rate, energy is 
transferred from the coherent oscillations of a to the radiation produced by the axion decay. 
The radiation produced thanks to the decay of the axion will quickly dominate the expansion 
and the second radiation-dominated phase will take place. 

The Bardeen potential prior to decay is given by Eq. (4.53) while, after the decay, its 
evolution equation (2.29) reduces to 

1. 



<5>" + m& + 2{H 2 + H')<5> - -V 2 $ = 0, v>Vd, 



(4.57) 



and the corresponding exact solution can be expressed as [26] 

$k(v) — — Bi{k){urj cos urj — sincj^) + B 2 (k) (cur/ sin cur/ + coscjt^J , i] > r^, (4.58) 

where u = k/y/3. In the sudden approximation, the two (dimensional) arbitrary constants 
Bi(k) and B 2 (k) can be uniquely fixed by matching, at the decay time rjd, the solutions 
(4.53) and (4.58) together with their first derivatives. The terms containing (kr] d ) are small 
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and negligible for modes that are outside the horizon at the time of the axion decay, i.e. for 
the ones relevant to the physics of the observed CMBR anisotropies. Furthermore, terms 
proportional to inverse powers of m/H& ~ Mp/m 2 are also small and can be consistently 
neglected. Hence, up to subleading terms, the final value of the Bardeen potential can be 
written as 



$ fc (7j) = $ (fc) 



2 (■( is ( — 



COS LOT] sin CUT] 

(urrj) 2 (urj) 3 



(4.59) 



where (3 = mr)da(r)d) ~ m/H& ~ Mp/m 2 . The /9-dependent prefactor is a consequence of 
the approximation of sudden decay where the axion field is assumed to decay at a specific 
time 77d- This sudden approximation also neglects the possible (exponential) damping of the 
oscillations in $ fc arising in Eq. (4.53). 

It will now be shown that the /^-dependent prefactor is an artefact of the sudden ap- 
proximation. In a realistic model of decay, in fact, the energy- momentum tensors of the 
radiation fluid and of the axion will not be separately conserved, because of their relative 
coupling induced by the friction term T(a'/a) 2 , which leads, in cosmic time, to the general- 
ized conservation equations: 

p a + (3H + r)(p a + Pa ) = 0, 

p r + 4#p r -r(p CT +p CT ) =0. (4.60) 

The fluctuations Xk will experience a similar damping, 

X " + (2H + Ya) X ' - V 2 X + -^a 2 X ~ V$' + 2^a 2 $ = 0, (4.61) 

while the Q k evolution will still be described by Eq. (2.29). This treatment of the damping 
of the fluctuations was suggested in [8] (see also [41]). The effect of T is, primarily, to induce 
a damping in the oscillations of the background and of the axion fluctuations according to 
Eqs. (4.60) and (4.61). Moreover, the (damped) fluctuations of the axionic field will also 
influence the dynamics of according to Eq. (2.29). The time-dependent oscillations of 
Eq. (4.53) (occurring in the absence of friction) will then be further suppressed if T ^ 
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(more details will be given in the following section). As a consequence, the /^-dependent 
correction tends to disappear from Eq. (4.59), leading to the final result 



Mv) = 3$o(fc) 



sin cur? cos ui] 



(4.62) 



(our/) 3 (ur)) 2 

In the equation for $ fc the effect of the finite duration T _1 is then equivalent to averaging 
over the decay time. At the end of the following section, numerical examples of the decay 
will be discussed in detail. 



V. BACKGROUND AND PERTURBATION EQUATIONS FOR a, > 1 

If cTj > 1, the epoch of axion domination precedes the oscillation epoch. The previous 
solutions for the slow-roll regime are still valid, and the axion starts dominating when 

p r a 2 = QH 2 ~ Va 2 , (5.1) 

i.e., for a quadratic potential, when 

(12)V4 

T = T a ~ 1-^. 5.2 

If <7i ~ 1, then 

r m ~ T a , (5.3) 

the axion oscillates almost immediately after becoming dominant, and the amplitude of the 
Bardcen potential at the onset of the oscillatory phase is obtained from Eq. (4.20) as 

Mv™) * -^Xi(k)a^ = ~ -I Xi (fc). (5.4) 

If the initial value a x is larger than 1, but not too large, then Eqs. (5.3) and (4.20) are still 
valid, but Eq. (5.4) is to be multiplied by the factor of, arising from a short period of axion 
dominance toward the end of the slow-roll evolution (see below). This effect, for moderate 
values of o"i, is illustrated in Fig 6 where, for the given parameters of the plot, the final 
amplitude of the Bardeen potential is estimated as 
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FIG. 6. The full curves are the results of numerical integrations for the case a x > 1. The dashed 
curves correspond to the approximated results of Eqs. (4.20) and (5.5). 



e 6 = 0.143 



(5.5) 



still in good agreement with the approximate value 1/7 of Eq. (5.4). 

If <J[ ^> 1, then r m ^> r a , and a phase of inflationary expansion dominated by the 
axion potential will take place between r a and r m . During this phase the axion slowly 
rolls, the radiation energy-density is quickly diluted as p T ~ a~ 4 , and the time evolution of 
the fluctuations is correspondingly modified. The Hamiltonian and momentum constraints 
(2.20) and (2.21) can now be combined to give 



4V 2 $ fc = 3Ha' Xk - $ k a' 2 + a' X ' k + V a a 2 Xk , 



(5.6) 



and the speed of sound of Eq. (2.33) becomes 



2a 2 



1 + 



(5.7) 



from which, using Eq. (2.40), 

2v; 



5p 



nad 



3Ha> 



$ k a' 2 + a' X ' k ~ 3Ha' Xk - a 2 V aXk 



(5.8) 



The combination of Eqs. (5.6) and (5.8) leads to 
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= ^V 2 $ k . (5.9) 

Hence, from Eq. (2.35), we get that (' k ~ at large scales. 

According to its definition, on the other hand, the constancy of (k implies, in cosmic 
time, that 

+ , 1 



+ (5.10) 



is also a constant, where 



ai = -H/H 2 . (5.11) 
It follows that, during inflation, we can parametrize the evolution of to lowest order, as 

$ fc = A(H/H 2 ), (5.12) 



where A is a constant controlled by the value of the Bardeen potential at the beginning 
of inflation. Assuming a quadratic potential for r < T a we have, from Eq. (5.4), A ~ 
— (l/7)(x\{k) / <J\)(H 2 / H) Ta . By using the dynamics of slow-roll inflation, 

H 2 * ^a 2 , (5.13) 

(5.14) 



m 



3# ' 
we can deduce that 

if 4 . . 

TP - (5.15) 

Using Eq. (5.15) we finally obtain the Bardeen potential at the onset of the phase of a- 
dominated oscillations: 

$ fe (r m ) ~<7iXi(A;), (5.16) 

where we used the expression for A given above and the fact that cr(r CT ) = <7i and cr(r m ) ~ 
0(1). For r > r m the axion eventually oscillates, and the subsequent evolution has the same 
features as already discussed in the previous section, for the case a x < 1. 
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In the following, the dynamics of the decay will be investigated numerically, and Eq. 
(2.29) will be solved together with Eqs. (4.60) and (4.61). In order to illustrate the results 
let us recall that, in the absence of friction (r = in Eqs. (4.61)), the evolution of Xk 
and during the axion-dominated oscillations, is given by Eqs. (4.53) and (4.54). In 
particular, 

<Z> k ( v ) ~ $ (k) + 5$ k (v), (5-17) 

where S^ k (rj) is an oscillating function* decaying as rj~ 3 ~ t' 1 . The frequency of oscillation 
of 5Q k is controlled by the axion mass. In analogy with Eq. (5.17) we can also define 5\k 
which, for T = 0, corresponds to the oscillating function appearing in Eq. (4.54). The 
evolution of t5Q k and of 5xk, for T = 0, is represented by the full bold curves of Figs. 7 
and 8. Notice that 5§ k and 5\k oscillate very fast and that, for our illustrative purpose, 
we have plotted their amplitudes calculated as the average of the semi- difference between 
the maximum and the minimum of each oscillation, and the semi-difference between the 
successive maximum and the same minimum. 

If the oscillations of are only suppressed by a power-law function of time, we have 
seen that there are mass- dependent terms that appear in the amplitude of the Bardeen 
potential after the decay. The integration of Eqs. (2.29) and of (4.60) and (4.61) shows 
however that, with the inclusion of the appropriate friction terms (due to the decay) into the 
energy-momentum conservation equations, the oscillations in 5§ k and 5\k are exponentially 
suppressed and, in such a case, no mass-dependent correction is left in the amplitude of 
the Bardeen potential (the asymptotic, constant value of is however unaffected by such 
a damping mechanism). The damping of the oscillations, on a time scale of order is 
illustrated by the dashed curves of Figs. 7 and 8. 

If we compare the sudden-decay approximation, discussed in the previous section, with 



*In order to avoid confusion we note that 5§ k and, in the following, 5xk, are not the power spectra 
of $ and x- 
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FIG. 7. Time evolution of the amplitude of the 5<&k oscillations (multiplied by t in cosmic 
time), with and without the damping term due to the axion decay. 

the numerical results of Figs. 7 and 8, we see that the finite duration of the decay process 
can be physically represented as a dynamical average to zero of the oscillatory terms in the 
evolution of In view of these results, when matching ^ to the post-decay phase, we 
should take into account the fact that all the derivatives of are exponentially suppressed 
with respect to $k/td, and thus can be safely neglected. This leads to the result reported in 
Eq. (4.62). 



VI. LARGE-SCALE ADIABATIC FLUCTUATIONS 

In order to discuss the direct impact of our results on the possible generation of the 
observed CMBR anisotropies, the evolution of the large-scale metric fluctuations should be 
followed down to the matter-dominated phase, for all times rj > i] cq . In particular, the phase 
and the amplitude of the Bardeen potential prior to i] cq will fix the initial conditions for the 
subsequent evolution of the inhomogeneities, and will be crucial to determine whether they 
are of adiabatic or isocurvature nature. 

We recall that, after the axion decay, the amplitude of the Bardeen potential has been 
computed as 
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FIG. 8. Time evolution of the amplitude of the 8\k oscillations, with and without the damping 
term due to the axion decay. 



* k {rj) = 3$ (fc) 



sin tor] cos cur] 



V < Veq, 



(6.1) 



(cur}) 3 [ur]) 2 J ' 

where, as in the previous section, uo = k/y/3. For rj > r/ eq , matter domination sets in, the 
background satisfies TH! + TC 2 = 0, so that the evolution of the Bardeen potential (outside 
the horizon) is described by 



+ m& k = o, n = -, 



(6.2) 



whose solution can be written as 



B(k) 

$ k (r)) = A(k) + T, > Tj cq . 



(6.3) 



Imposing the continuity of the solutions (6.1) and (6.3) and of their first derivatives at 
V = ^eq one obtains: 



A{k) = 



B{k) 



3$ (fc) 
3%(k)r), 



[-2x cq cosx cq + (x 2 + 2) sinxeq], 



cq 



5x^ q 



[3;c e q cosx e q 4~ (x eq 3) sinx e q], 



(6.4) 
(6.5) 



where x eq = ojr] cq = kr] cq /\^3. For scales that are outside the horizon prior to decoupling, 
x cq <C 1, and Eq. (6.3) becomes 
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Mv)=Mk)[l+ { -^(^) 5 ], V>Vo q - (6.6) 

For i] > r/ eq the decaying mode is highly suppressed, and we are then in the situation of 
constant Bardeen potential right after equality with an amplitude &o(k), which (recalling 
the previous results (4.43), (5.4), (5.16)) is completely determined by the axion spectrum 
and by the initial conditions of the axion background. More precisely, the final amplitude 
can be parametrized as follows 

= MVd) = -/(^)Xi(A:), (6.7) 

where 

/(<7i) = Ci<7i + — - c 3 , (6.8) 

and 

ci ~ 0.13, c 2 ~ 0.25, c 3 ~ 0.01. (6.9) 

The above coefficients q have been obtained by integrating numerically the evolution equa- 
tions of the background and of the fluctuations for different values of <7i (both larger and 
smaller than 1). Then, following the hint of the analytical results obtained by solving the 
evolution piece-wise, the final value of ®k(v) nas been fitted with Eq. (6.8), and the values 
reported in Eq. (6.9) have been determined. 

The value (6.8) of the Bardeen potential provides the initial condition for the subse- 
quent hydrodynamical evolution. Such evolution will allow us to determine, in turn, the 
precise value of the temperature fluctuations through the Sachs- Wolfe effect. In particular, 
the modes that are outside the horizon for i] cq < r] < r^cc will determine the large-scale 
temperature fluctuations relevant to the COBE observations. 

By perturbing the corresponding conservation equations on a matter-dominated back- 
ground, we obtain: 

K ~ ^VV - 4$' = 0, (6.10) 
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<-^r-$ = 0, (6.11) 

S' m -V 2 v m -3& = 0, (6.12) 
v' m + Hv m -$ = 0, (6.13) 

where 8 m = Sp m /p ni and v m , following the notation of the previous sections, are the gauge- 
invariant density contrast and velocity potential of the matter fluctuations. Also, in the 
above equations, 

2 24 

?{ = -, PraCL 2 = —. (6.14) 
7] I] 2 

As already stressed at the beginning of this section, <£> is constant during the matter- 
dominated phase. Using this property we can now work out the specific relations between 
the different fluid variables, for modes that are outside the horizon right after equality, so 
as to explicitly check the adiabaticity of the fluid perturbations. 

The system of Eqs. (6.10)-(6.13) can be easily solved by going to Fourier space. For v m 
we have 

kv m (k)~^-$ (k), to7«l. (6.15) 

Since Vi> m (evaluated outside the horizon) contributes directly to the Sachs- Wolfe effect, it 
is important to notice that this term is subleading with respect to the other contributions 

— * 

arising in the case of adiabatic fluctuations. We will indeed show that, unlike Vf m , which 
is suppressed, the contrast 5 r is instead constant outside the horizon, and proportional to 

^(AO- 
Inserting v T from Eq. (6.10) into Eq. (6.11) we get a decoupled equation for 8 T , namely, 

5>: + j5 T = ~k 2 %(k)- (6.16) 

The general solution is 

5 r (krj) = A 1 coslut] + B 1 sinurj + 4& (k)[cosLur] — 1], (6-17) 
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and the constants A\ and B\ can be determined by consistency with the other equations 
and with the Hamiltonian constraint (2.20) written in the case of a matter-radiation fluid. 
The final result is: 

4 r i 
S T (k,rj) = -$ (fc) coslut] - 3 (6.18) 
3 

kv T (k,rj) = sin uitj, (6.19) 

V3 

5 m (k,r,) = -2* (k)-^p-(kn) 2 , (6.20) 
fcw m (fc,i 7 ) = M$ ( fc ). (6.21) 

Notice that, outside the horizon, kv m = kv T as required by local thermodynamical equilib- 
rium. Furthermore, for krj <C 1, the velocities of the two fluids are proportional to (krj). 
When the modes are outside the horizon, Eqs. (6.18) and (6.20) imply that the density 
contrasts 5 r and 5 m are both constant and proportional according to 

5 r ~ (4/3)<J m . (6.22) 

This result has a simple physical interpretation, and implies the adiabaticity of the fluid 
perturbations. The entropy per matter particle is indeed proportional to S — T 3 /n m , where 
n m is the number density of matter particles and T is the radiation temperature. The 
associated entropy fluctuation, 5S, satisfies 

$S 3 „ . 

-j = fr ~ S m , (6.23) 

where we used the fact that p r ~ T 4 and that p m = mn m , where m is the typical mass of 
the particles in the matter fluid. Equation (6.22) thus implies SS/S = 0, in agreement with 
the adiabaticity of the fluctuations. 



A. Sachs- Wolfe effect and COBE scales 



The fluctuations of the Bardeen potential and of the radiation density contrast are sources 
of a slight temperature difference between photons coming from different sky directions. 
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This is the essence of the Sachs- Wolfe effect [42] . In terms of the gauge-invariant variables 
introduced in the present analysis, the various contributions to the Sachs- Wolfe effect, along 
the n direction, can be written as [26,29] 

AT [5 - i f r ' dcc 

—(n,r) ,x )=j + n-Vv b + <S>{r) dec ,x{r) dec ))- ($' + V')(ri, x(r)))dr), (6.24) 

where i] is the present time, and x{rj) = x — n(r] — ri ) is the unperturbed photon position at 
the time rj for an observer in xq. The term v b is the peculiar velocity of the baryonic matter 
component. We are preliminarily interested in the effects of scales still outside the horizon at 
the time of the matter-radiation equality, which are the scales relevant to the observations of 
the COBE-DMR experiment [14,43]. In order to correctly take into account the constraints 
imposed by the COBE normalization on the spectral amplitude of the Bardeen potential, let 
us compare the relative weight of the different terms appearing in the Sachs- Wolfe formula 
(6.24). 

From Eq. (6.19) we can see that, for our adiabatic initial conditions, the fluctuation in 
the matter velocity potential is subleading for superhorizon scales, suppressed by the term 
ki] <C 1 with respect to the constant values of 5 r and Furthermore, since $' ~ and 
\& = <E», the integrated Sachs- Wolfe effect can also be neglected. By inserting Eq. (6.18) 
into Eq. (6.24) we thus obtain the usual result for adiabatic fluctuations, namely 

AT 1 

— (n, rj , x ) = -$(?7dec, x(rj dec )), (6.25) 

to be used for the comparison of our theoretical predictions with the COBE normalization. 

On the other hand, by taking the Legendre transform at the present time 770, the temper- 
ature fluctuations of Eq. (6.24) can be generally expanded into spherical harmonic functions, 

Yim, as 

AT 

— (x , n, 770) = ^2 a em(x )Y em (n), (6.26) 

£,m 

where the coefficients ap m define the angular power spectrum Cp by 

dim ■ de'm'} = die'fimm'Ce, (6.27) 
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and determine the two-point correlation function of the temperature fluctuations, namely 
/ ST 6T \ 

(jr(n) — (n)\ = ^ ^ae m a* eim ^Y em (n)Y e 1 ml (n') 

' ' (n-3'=costf) tt'mm' 

= y J2( 2i + l)C^(cos^). (6.28) 

These coefficients Ce, in turn, are related through Eq. (6.25), to the power spectrum of 

&o(k), and for 2 < £ <C 100 they can be expressed as [44] 

9 f°° r\h i \ 
C '-^J ^(\Mk)\ 2 )k 3 mVo-mec)}. (6.29) 

As already stressed, the spectrum of the Bardeen potential is fully determined, in our 
context, by the initial spectrum of axionic fluctuations amplified by the pre-big bang dy- 
namics. A self-contained derivation of such a spectrum, including the mass contribution, is 
presented in Appendix A. Consider first the case of minimal pre-big bang models, whose 
related spectrum is reported in Eq. (A. 14). The spectrum of curvature perturbations will 
then be, at large scales, 

k*\%(k)\ 2 = f\a^\ Xk \ 2 = f 2 ^^) \±J , k < hi, (6.30) 

where k\ is the maximal amplified comoving frequency, i.e., in our conventions, the frequency 
at which only one axion per cell of phase-space is produced. From hi, a typical curvature 
scale Hi (which can be, at most, of the order of the string mass) can be obtained: 

H 1 = — < M s . (6.31) 
ai 

The particular value of the scale Hi may be regarded as a phenomenological parameter of the 
chosen model of pre-big bang evolution. Even assuming, according to the standard lore [45], 
that M s ~ 10 _1 Mp, still the exact relation of H 1 to M s depends on the detailed dynamics of 
a highly curved and strongly coupled background. The approach of the present investigation 
has been to include all the theoretical indetermination into Hi, trying to have a reasonable 
control of all the other numerical factors associated with the post-big bang evolution. In Eq. 
(6.30) the particular value of n depends upon the specific model of pre-big bang evolution 
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[15,21]. In the case of a ten-dimensional model with an isotropic six-dimensional internal 
space, the line element can be written as 

ds 2 = dt 2 - a 2 (t) lij dx i dx j - b 2 {t) lab dy a dy\ (6.32) 

where i,j run over the three external space-like dimensions and a, b run over the six internal 
dimensions. Defining as 

r = (6.33) 
2V 6 V 3 

the relative rate of variation of the external V 3 = a 3 and internal V 6 = b 6 volumes, the 
spectral index n can be expressed as [21] 



4 + 6r 2 - 2V3 + 6r 2 , r 

n = . 6.34 

1 + 3r 2 y ' 

The case of flat spectrum (i.e. n — 1) corresponds to the case r = ±1. If internal dimensions 

are static (i.e. r = 0), then n = 4 — 2\^3 ~ 0.53. Blue spectra are allowed when the rate 

of variation of the external volume is much smaller than the internal one. The maximal 

n achievable in this case is n = 2, corresponding to the case of static external manifold 

(r — > oo). 

Bearing in mind Eqs. (6.31) and (6.34), we can use Eq. (6.30) and perform the integral 
of Eq. (6.29). For — 3 < n < 3 the integral appearing in Eq. (6.29) can be done analytically 
[44] and the result is 



c ( SW) _2 n 2( /^x 2 ro y-i r(3-n)r(f-| + f) 

Here c^o — 10~ 18 Hz and ui(to) = Hiai/a are, respectively, the proper frequencies corre- 
sponding to the present horizon scale and to the present value of the cut-off scale ki. 

wi(fo) = (H ± iJ cq ) 1 / 2 (I) 1/6 ar 2/ ^e" q 1 , * < 1, (6-36) 

u^to) = (H, H CCi f 2 (-) l '\] l2 Z- l z-^ a, > 1. (6.37) 



(we have rescaled taking into account the kinematics of the various cosmological phases 
from t\ down to to). The factor Z a = (a osc /a a ) denotes the amplification of the scale factor 
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during the phase of axion-dominated, slow-roll inflation, for the case a\ > 1. Notice that 
a>i(t ) depends on the mass, on the initial amplitude of the axion background and on the 
axion decay rate. If the axion decays at a typical scale fixed by Eq. (3.4), Eqs. (6.36) and 
(6.37) lead to 

(we have used ff — 10~ 6 H eq ~ lCT 60 Mp). Hence, in spite of the fact that the initial axionic 
spectrum does not have any mass dependence, the mass appears again when computing the 
amplitude of the spectrum at the present horizon scale u . 

The amplitude of the Bardeen potential, on the other hand, is constrained by the COBE 
normalization of the quadrupole coefficient C2, which in our case is given by 

where 

„ 2 n T(3-n)r(^) 
a = — v 2 ' 6 41) 

n 72r 2 (^)r(^)' { } 

Using the experimental result [46] 

C 2 = (1.9 ±0.23) x 1(T 10 , (6.42) 

we are thus led to the bounds 

(5-n)/2 / \ -(n-l)/3 



jk) lO- 29 ^" 1 ^ 1-9 x 1(T 10 , OKI, (6.43) 

/ tt \ (5-n)/2 / \ -(n-l)/3 

^(aO^-Vf 1 -^ lO" 29 ^- 1 ) ^ 1.9 x 1(T 10 , a, > 1. (6.44) 



These constraints, imposed by the COBE normalization, will be discussed at the end of the 
present section, and combined with other theoretical constraints pertaining to the various 
models of background evolution. 
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B. Acoustic peak region 



In the previous discussion of the modes that are outside the horizon before decoupling, 
we have completely neglected the possible scattering of radiation with baryons. In fact, 
if we move to smaller angular scales (i.e. typically to t > 100), the main contribution 
to the CMBR temperature fluctuations comes from the oscillations of the various plasma 
quantities, the so-called Sakharov oscillations [47]. A correct approach to this problem is 
then to perturb consistently the Boltzmann equations for the different species of the plasma 
[48-50]. Furthermore it can be relevant to discuss the case of a smooth transition between 
radiation and matter dominated epochs. In such a context it becomes difficult to provide an 
analytical description of the system and, in order to compute the patterns of the acoustic 
oscillations, we will indeed present some numerical examples in the third part of the present 
section. 

It is however useful to emphasize that the phases of the Bardeen potential for the adia- 
batic mode of Eq. (6.1) determine not only the relative weight of the Sachs- Wolfe contribu- 
tions, but also the specific phase of the oscillatory patterns at small scales in the temperature 
fluctuations. For scales t ?t 100 the contribution to the temperature perturbations given in 
Eq. (6.24) is dominated by acoustic oscillations. This aspect can be appreciated by looking 
at Eqs. (6.18)-(6.21) in the limit krj > 1, where the peculiar velocity of baryonic matter 
does not oscillate. Instead, from Eq. (6.18), we find that the terms <5 r /4 and $, appearing 
in Eq. (6.24), combine to give a single term oscillating like a cosine: 

^(k,f] ,i] dcc ) ~ ^5 r {k,i] dcc ) + $ (fc) ~ cos ui] dec . (6.45) 

In this argument the interactions of baryons with the radiation fluid have been neglected. 
The dynamics of (AT/T)^ can be obtained from an exact Boltzmann equation with source 
term provided by Compton scattering coupled to the continuity and Euler equations for the 
fluid variables. Before recombination, Compton scattering is very rapid and therefore the 
Boltzmann, Euler and continuity equations for the photon-baryon system can be expanded 
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in powers of the Compton scattering time [49,50]. Within this approximation the baryon 
velocity field is damped and (AT/T) k oscillates as a cosine for adiabatic initial conditions. 
In the approximation of [49,50], the oscillations in (AT/T)^ have an amplitude proportional 
to (1 + i?)~ 1 / 4 where R(rj) = 3pb/(4p r )- This result simply tells that the baryonic content of 
the plasma determines the height of the first peak. Notice that this is in sharp contrast with 
what happens in the case of light axions [17,18], where the Bardeen potential is quadratic in 
the axion fluctuations, and the initial conditions for the hydrodynamical evolution are of the 
isocurvature type. This implies, in particular, that the oscillatory patterns of the CMBR 
anisotropies will be shifted by n/2 if compared with the case discussed in the present paper. 

C. Constraints on pre-big bang models 

In this subsection we will discuss the bounds imposed by the COBE normalization, 
together with other constraints following from the evolution of the background geometry. 
Let us start with the axion spectrum of minimal pre-big bang models, Eq. (6.30). In such 
a case, and for a flat Harrison-Zeldovich spectrum (i.e. n = 1), the COBE normalization is 
inconsistent with a cut-off Hi at the standard value M s ~ 10 _1 M P of the string mass scale 
[45]. By using n — 1, and taking for U; the value minimizing /(o"i), 

of* = ^ ~ 1.38, f{vf n ) ^ 0.34, (6.46) 

we have indeed, from Eqs. (6.40)-(6.42), 

H x ~ 5.2 x 1(T 4 M P . (6.47) 

However, the precise value of Hi is one of the main uncertainties of pre-big bang models. As 
we shall see in a moment, the value Hi = M s (or Hi = Mqut) m ay become consistent with 
the COBE normalization for non-flat (blue) spectra, and even for a strictly flat spectrum in 
the case of non-minimal implementations of the pre-big bang scenarios. 

Let us first recall the various constraints to be imposed on the spectrum. The condition 
(6.43) is to be combined with the constraint (3.7), the condition (6.44) with the constraint 
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(3.12), which are required for the consistency of the corresponding classes of background 
evolution. Both conditions are to be intersected with the experimentally allowed range of 
the spectral index. We will use (as a reference value) the generous upper bound [1], n < 1.4. 
Also, for our illustrative purpose, we will take the maximally extended range of allowed 
values of the axion mass, satisfying the nucleosynthesis constraint m > 10 TeV. 

We will assume, finally, that in the case <7i > 1 the axion-driven inflation is short enough, 
to avoid a possible contribution to Ce arising from the metric fluctuations directly amplified 
from the vacuum during such a phase of axionic inflation. This requires that the smallest 
amplified frequency mode u a , crossing the horizon at the beginning of inflation, at decou- 
pling be still larger than the Hubble horizon at the corresponding epoch. This imposes the 
condition u a (t ) = H a (a a /a ) > u dcc (t ) = H dcc (a dcc /a ) , namely 



to be added to the constraint (3.12) for o x > 1. 

The allowed region in the plane {log<7i, log(m/M P )} is illustrated in Fig. 9 for Hi = 
lCT 2 Mp, using for the inflation factor the parametrization Z a = exp((of — l)/8). Along 
the thin full curves the parameters satisfy the COBE normalization, for fixed values of n, 
ranging from 1.1 to 1.4 (the condition (6.48), in this case, is always automatically satisfied). 
A growing ("blue") spectrum is thus allowed even if Hi ~ M s , for a wide range of axion 
masses, and for a (narrower) range of values of (T;. In particular, for the case Hi = 10~ 2 M P , 
we find 1 > o"i > 10~ 4 , for o x < 1. For o"; > 1 the results are complementary for the spectral 
index, but there are much more stringent bounds on <7i, because the inflationary red-shift 
factor Z a grows exponentially with of. As a consequence, the allowed region for a x > 1 is 
distorted and compressed, as illustrated in Fig. 9. 

The allowed region may be further extended if the inflation scale Hi is lowered (see 
for instance [51]), and a flat (n — 1) or almost flat spectrum may become possible if 
c 2 aiHi < 10~ 5 M P (Ti, for a x < 1, and if ciCtiHi < l(T 5 Mp/(Ti, for o"; > 1 (see Eqs. (6.43), 
(6.44)). The corresponding allowed values of H x and a x are illustrated in Fig. 10 for 




(6.48) 
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log C; 

FIG. 9. Allowed values of a\ and m (in Planck units) according to Eqs. (6.43), (6.44), with 
ill = 10 _2 M P . The allowed region (within the thick lines) is bounded by the condition n < 1.4 
(left and right bold lines), by the nuclesosynthesis constraint m > 10 TeV (lower bold line), by the 
condition (3.7) (upper left bold line) and (3.12) (upper right bold line). 

m = 10 _9 Mp, and for three different values of n around 1. 

However, a flat spectrum may be allowed even keeping pre-big bang inflation at the 
string scale (Hi ~ M s ), provided we consider a non-minimal pre-big bang scenario. In that 
context, in fact, the high-frequency branch of the axionic spectrum may be modified, getting 
steeper enough to match the string-scale normalization at the end-point of the spectrum, 
while the low-frequency branch remains flat (or quasi-flat, see Appendix), to agree with 
large-scale observations. Examples of realistic pre-big bang backgrounds producing such an 
axion spectrum have been presented already in [19]. 

A non-minimal spectrum can be parametrized by the Bogoliubov coefficients (which will 
be given in Eq. (A. 16)), in terms of a generic break-scale k s and of the high-frequency 
slope parameter 5. In that case, for a long and/or steep enough high-frequency branch 
of the spectrum, the large-scale amplitude may be suppressed sufficiently to allow flat (or 
even red) spectra at the COBE scale. In fact, for the non-minimal spectrum (A. 16), the 
normalization condition (6.40) becomes 
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FIG. 10. Allowed values of H\ as a function of crj for different values of the spectral index and 
for m = l(T 9 Mp. 

Flat or red spectra (n < 1) are thus possible even for #! > KT 2 M P , provided 

) < 10- 6 . (6.50) 
^>\) 

In order to illustrate this possibility we will choose a specific model of background by 
identifying k s with the equilibrium scale k eq , in such a way that n corresponds to the spectral 
index of all scales relevant to the CMBR anisotropics, while n + 5 provides the average 
spectral index for all other scales, up to k\. We will also assume for the axion background 
the "natural" initial value a\ — 1, so that 

w i _ . ., ^27 ( ( m \ 1/3 



— 10 7T" ) hnr ■ (6-51) 
a; s cu eq VMp/ VM P / 



The COBE normalization can then be written explicitly as 

By using the experimental value of Ci given in Eq. (6.42) we can now obtain a relation 
between the high-frequency slope parameter 5 and the spectral index n at the COBE scale, 
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FIG. 11. Relation between 5 and n for different values of H\ (in Planck units), for m = 10 _9 Mp, 
and for a\ = 1. The vertical dashed lines denote the experimentally allowed range 0.87 < n < 1.06. 

for any given value of Hi and m. In Fig. 11 we illustrate such a relation for different 
(realistic) values of Hi, and for a typical axion mass m = 10 _9 Mp. It should be stressed 
that, for n ~ 1, and a x of order 1 (i.e. near the minimum of /(o"i)), the curves at constant Hi 
are almost insensitive to the values of m, and remain stable even if we change m by various 
orders of magnitude, as illustrated in Fig. 12. 

We have also reported, in Fig. 11, the (present) most stringent bounds on n, obtained 
by a recent analysis of the CMBR anisotropies and large-scale structures [55,56], i.e. 0.87 < 
n < 1.06. They are all compatible with Hi ~ M B , provided we allow for a small break of the 
minimal spectrum, with 5 ~ 0.2-0.3. On the other hand, as already stressed, no break at 
all is needed (i.e. 5 = 0) if, for some dynamical mechanism (see for instance [51]) the string 
mass is lowered down to the GUT scale, i.e. Hi ~ 10 _3 Mp. 

In Fig. 12 we have plotted the same curves of Fig. 11 for two, very different values 
of the axion mass, 10~ 9 Mp (bold curves) and 10~ 14 Mp (thin dashed curves). As clearly 
illustrated by the figure, the dependence on the mass is very mild, and it becomes practically 
inappreciable (for the given range of parameters) when Hi approaches Mqut- 

Having discussed the constraints imposed by the COBE normalization we can present 
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FIG. 12. Stability of the curves of Fig. 10 for two different choices of the axion mass, 10~ 9 Mp 
(bold curves) and 10~ 14 Mp (thin dashed curves). 

now the plots of the angular coefficients Ci for the scalar spectrum, in the case of a spatially 
flat background and for few relevant choices of the spectral index. In Figs. 13 and 14 the 
(scalar) angular power spectrum defined in Eq. (6.27) is reported for a flat, slightly red 
and slightly blue spectrum. In order to obtain the results of Figs. 13 and 14 we used the 
latest release of CMBFAST, relaxing the strict COBE normalization at I = 10, in favour of 
a better general agreement of the overall fit. 

In Figs. 13 and 14 we used the following values of the cosmological parameters, selected 
according to fits of CMBR anisotropy experiments [56]: h = 0.65, — 0.7, and h^fli, = 
0.02. The selected value of h^Qb is rather robust, even if no final consensus has been reached 
on the second significant figure beyond 0.02. We have also assumed the simplest scenario 
for the late-time cosmological evolution, with no significant effects of reionization. 

In Fig. 13 the Ci are plotted on a linear scale, whereas in Fig. 14 we present the 
same plot with a semi-logarithmic scale, in such a way that the region relevant to the COBE 
observations is less compressed. We recall, finally, that the flat, red and blue spectral indices 
may correspond to particular combinations of the parameters Hi,m,S and a\, chosen in such 
a way as to satisfy the COBE normalization, Eq. (6.49). 
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FIG. 13. The spectrum of Ci is illustrated for a fiducial set of parameters (ho = 0.65, 
Q h = 0.04733, n A = 0.7, n m = 0.25267) and for flat (full line, n = 1), slightly red (dashed 
line, n = 0.9) and slightly blue (dotted line, n = 1.02) spectral indices. 

The data points reported in Figs. 13 and 14 are those from COBE [1,2], BOOMERANG 
[4], DASI [5], MAXIMA [6] and ARCHEOPS [7]. Notice that the data reported in [7] fill 
the "gap" between the last COBE points and the points of [4-6]. Therefore, one could 
think of normalizing the spectra not to COBE but directly to ARCHEOPS. In spite of this 
possibility, the forthcoming MAP data will give even more accurate determination of the Ce 
spectra. It will then be interesting to use these data in order to make more consistent and 
accurate determinations of the pre-big bang parameter space. 



VII. CONCLUDING REMARKS 

In the present paper the possible conversion of isocurvature, primordial axionic fluctu- 
ations into adiabatic, large-scale metric perturbations has been discussed in the context 
of the pre-big bang scenario. Depending upon the specific relaxation of the axionic back- 
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FIG. 14. The same plot as in Fig. 13, but on a semi-logarithmic scale. 

ground toward the minimum of the potential, a constant (and large enough) mode in the 
Bardeen potential can be generated, for scales that are still outside the horizon right after 
matter-radiation equality. 

After analysing the dynamics of the background and of its fluctuations, the final ampli- 
tude and spectrum of the Bardeen potential has been related to the initial axion spectrum 
directly arising from the vacuum fluctuations amplified during the pre-big bang epoch. Our 
goal has been to include, with reasonable accuracy, the details of the post-big bang evolu- 
tion, in such a way that the pre-big bang parameters could be directly constrained by the 
COBE normalization, and by the analysis of the Doppler-peak structure. All the theoretical 
uncertainty reflects in our lack of knowledge of H\ which determines the end point of the 
primordal axion spectrum. 

The main conclusion of this work is that a phenomenologically appealing spectrum of 
adiabatic scalar perturbations can naturally emerge from the simplest pre-big bang scenario 
through a conversion of the initial isocurvature perturbations of the Kalb-Ramond axion. 
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Since, at the large scales tested by CMBR experiments, the above conversion preserves the 
scale-dependence of the original spectrum, it is important for the latter to be quasi-scale- 
invariant at large scales. This can be achieved, for instance, if the very early stages of pre-big 
bang cosmology at weak coupling involve a symmetric evolution of all 9 spatial dimensions 
(modulo T-duality). Since the constant mode of the curvature fluctuations leads to adiabatic 
initial conditions for the fluid evolution after matter-radiation equality, the location of the 
Doppler peak is correctly reproduced. 

On the other hand, the absolute normalization of fluctuations at large scales (say those 
relevant for COBE) depend on several details of the model. Indeed, the axion spectrum 
is naturally normalized at its end-point, given by our parameter Hi. If one takes, naively, 
Hi ~ M s ~ 10 17 GeV and assumes a flat spectrum one finds values of AT/T that are a 
couple of orders of magnitude too large when compared with COBE's data. However, one 
can think of many (individual or combined) effects that can bring down our normalization 
to agree with the data, e.g. 

• A slight (blue) tilt to the spectrum; 

• A blue spectrum just at high frequency (i.e. for scales that exit late, during the 
strongly coupled regime); 

• A lower Hi/M s ratio; 

• A lower M s /Mp ratio. 

In the near future we hope to extend the present discussion to forthcoming CMBR 
anisotropy data at even smaller angular scales. It would be interesting to see if a combined 
analysis of the experimental data may give further useful hints on the parameter space of 
the scenario explored in the present investigation. 
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APPENDIX A: AXIONIC SPECTRA 



During the pre-big bang phase the quantum mechanical fluctuations of the axionic field 
will be amplified from the initial vacuum state. The obtained spectrum provides the initial 
condition for the evolution of the axion fluctuations in the post-big bang phase. At very 
large-scales, such a spectrum will not depend so much upon the details of the pre-big bang 
evolution. At smaller scales, however, it can be strongly affected by specific dynamics of the 
strong coupling and high-curvature regime. In spite of the fact that the spectral slope at large 
scales is not affected by high energy corrections, the large scale amplitude is affected and, in 
particular, a steeper slope at small scales has important consequences for the normalization 
of the low-frequency branch of the spectrum. In this appendix we will consider, separately, 
the axion spectrum obtained in the case of minimal and non-minimal pre-big bang models. 



1. Minimal pre-big bang models 

The linearized evolution of massive axion inhomogeneities Xk, neglecting their coupling 

to scalar metric perturbations, in a spatially flat cosmological background, is described in 
general by the equation 

\k 2 + m 2 a 2 -— U k = 0, (A.l) 
L z J 

where 

z = a e^ 2 , ip k = zxk- (A. 2) 

In the pre-big bang phase (77 < 771) the axion is massless. In the post-big bang, radiation- 
dominated phase, taking place for 77 > rji, the gauge coupling freezes {ip = const) and the 
axion acquires a mass. The produced axion spectrum, in principle, has a relativistic and 
a non-relativistic branch: this is because, in the radiation era, the proper momentum is 
red-shifted with respect to the rest mass, and the whole spectrum, mode by mode, tends to 
become non-relativistic. The spectral slope of the relativistic and non-relativistic branches 
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of the spectrum are in general different. However, if the axion modes, as in the present case, 
become non-relativistic when they are still outside the horizon, the solution is then exactly 
the same as in the relativistic limit. 

Consider first the relativistic branch of the spectrum. For rj < rji the solution of Eq. 
(A.l) can be expressed in terms of the second-kind Hankel functions [40] as: 

Mv)=V 1/2 H^(kr ] ), (A.3) 

where \i depends on the parameters controlling the kinematics of the pre-big bang back- 
ground (a specific example will be given below, see Eqs. (A. 7) and (6.34)). In the radiation 
era, rj > rji, one has z" / z = 0, and the evolution equation of ipk acquires a massive correction: 

< + (k 2 + m 2 a 2 ) 4> k = 0, (A.4) 

Assuming that the axion mass is negligible at the transition epoch rji, the solution (A.3) 
can be matched to the plane-wave solution 

4k = ^= [c + {k)e-^ + c4k)e^] , (A.5) 

and the final result for Xk is 

xM = ^7rMkn\ (A.6) 

ay k 

where 

1 re — 5 

c W *(-)', (A.7) 

with n — 4 — 2\/jl\ . Note that the expression of the Bogoliubov coefficient c(k) and of the 
mean number of produced axions, = \c(k)\ 2 , contains different numerical factors of order 
1. At the same time the maximal amplified momentum k\ can be defined in different ways, 
all equivalent up to numerical factors. In the present analysis we will define the maximal 
scale k\ as the energy scale where one axion is produced per unit volume of phase-space. 
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Consider now the non-relativistic spectrum in the case when the mode becomes non- 
relativistic while it is still outside the horizon. Defining as k m the limiting comoving fre- 
quency of a mode that becomes non-relativistic (k m = ma m ) at the time it re-enters the 
horizon (k m = H m a m ), we find, in the radiation era [18,19], 

(in \ ^^ 2 
H~J ' (A ' 8) 

We are thus considering modes with k <C k m . In order to estimate the spectrum, in this 

limit, let us write Eq. (A. 4) in a form suitable for comparison with known results of parabolic 

cylinder equations: 

d 2 ^ k , /V_ 6 V fc = 0> x = n (2ayi 2 , -b = k 2 /2a, (A.9) 



dx 2 V 4 
where 

m 2 a 2 = a 2 r] 2 , a = mH\a\, (A. 10) 

and where a ~ 77 has been assumed. The corresponding general solution can be written as 

1> = A yi (b,x) + By 2 (b,x) , (A.ll) 

where y\ and y 2 are the even and odd parts of the parabolic cylinder functions [40]. The 
normalization to Eq. (A. 6) in the relativistic limit (i.e. x — > 0) gives A = and 

k n 1/2 



^ k ~ ° {k) m{b ' x) - (A ' 12) 

Outside the horizon, and for non-relativistic modes, k <C ma, we take (respectively) 

the limits — bx 2 <C 1 and -&<i 2 , the solution can be expanded as y% ~ x ~ r]\^2a, so that 
the mass disappears from the amplitude: 



\c(k)\k 1/2 

Xk ~ LU1 . (A.13) 

am 



The insertion of the spectrum (A. 7), using k\ = a\Hi, leads to the final result 



k*' 2 \ Xk \^ Hl (^y . (am) 
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2. Non-minimal pre-big bang evolution and spectral breaks 



Equation (A. 14) holds in the case of minimal pre-big bang models, where the dynamical 
evolution of the dilaton field is dictated by the solution of the low-energy equations of 
motion. However, when the dilaton enters the strong coupling regime, different types of 
scenarios may emerge. In particular, relation (A. 2) defining the form of the axion pump 
field, may change in the infinite bare string-coupling limit, as suggested by the arguments 
recently developed in [51]. In the framework of [51] the axion coupling function, as well as 
the other coupling functions pertaining to fields of different spin, may have a finite limit for 
infinite bare string coupling. Hence, toward the end of the pre-big bang phase (i.e. when 
strong coupling is presumably reached), 



where c z is a constant. Since the axionic pump field now depends only on the scale factor, 
it will naturally be steeper for small length scales. A complementary possibility, discussed 
in [19], is the presence of an intermediate high-energy phase, which precedes the standard 
radiation era, and which is still part of the accelerated pre-big bang regime, but in which the 
kinematics of the (usual) canonical pump field is significantly different from its low-energy 
behaviour. 

In all these cases the obtained spectra, at small scales, are possibly steeper than in the 
case of minimal pre-big bang models. In the simplest case the spectrum will have only one 
break, at a momentum scale that will be conventionally denoted by k s , and the Bogoliubov 
coefficients can be written in the form 



Here 5 > parametrizes the slope of the break at high frequency, while n is the usual spectral 
index appearing at large scales and computed on the basis of the perturbative evolution of 



z~a[c z + 0(e-^ 2 )], 



(A.15) 





(A.16) 
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the dilaton field. From Eq. (A. 16) it can be argued that the steeper and/or the longer the 
high-frequency branch of the spectrum, the larger the suppression at low-frequency scales. 
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